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A large class of classical dynamical systems with an external rapidly oscillating driving action 
is considered and the effective Hamiltonian-like equations for the mean motion are obtained. The 
respective Liouville equation for the distribution function of the mean coordinates and momenta is 
derived. PACS numbers 03.20.+i, 05.20.-y, 05.20.Gg 
The inverted pendulum, studied by P.L. Kapitza and recently revisited by S.-Y. Kim and B. Hu f|, is a very good 
example of systems where a fast driving external action imposes a complex dynamics to the system. The main result 
in Kapitza's theory is that a simple nonlinear oscillator under the action of a rapidly oscillating force is not conserving 
the mean motion on trajectories of a non perturbed analog. Instead of this, it behaves in a completely different 
manner, because an additional non small term proportional to the squared amplitude of external pulses appears in 
' the restoring force. This additional restoring force is responsible for the stabilization of the inverted pendulum, which 
, appears as a "miracle " and in no way could be predicted intuitively. It is therefore very tempting to look for such 
effects in extended many-body systems. 

In this Letter, we propose an approach to the treatment of a large class of such systems in a Hamiltonian-like 
formalism, which provides a straightforward transition to the statistical thermodynamics of these systems. 
Consider the following dynamical equations 
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where g and h are arbitrary functions of all canonical variables. Here Q 3> lo , where 2tt/lo is the characteristic 
time of the undriven system (with the Hamiltonian i?(q, p,t)). The series in the right-hand side of (1), (2) can be 
considered as Fourier expansions of arbitrary vector functions of t with characteristic time 2n/Q. 

We do not require the system being necessarily Hamiltonian, i.e. the equations (1), (2) are of the general type 
r= J~{T), where T (t) — (q, p) and T is a vector function of the phase space. The non-Hamiltonian nature of these 
t-H , equations means that the compressibility does not in general vanish: 
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0^ • The non-Hamiltonian dynamics is useful when considering open or driven systems, under the action of some external 
influence, such as heat bath or mechanical piston, etc. M. E. Tuckerman, C.J. Mundy and MX. Klein ||, jl| have 
recently derived a correct generalization of the Liouville equation to account for nonvanishing compressibility of phase 
space (3): 
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where the Jacobian J satisfies the differential equation 



Note that substituting (5) into (4) and taking into account the obvious relation ^ = T VrJ leads to the 

conservation condition df /dt = or 
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^+r-V r / = 0. (6) 

The backbone of our approach is the search of solutions of eqs. (1), (2) in the form 

qi = Qi+MXi! (?) 
Pi = Pi + p,pi, (8) 

where Q, P and Xi P are respectively the "slow " and "fast " parts, whose characteristic times are accordingly 
T ~ 2ir/u and r ~ 27r/il, and p = uj /fl <C 1. The expression (7) was successfully used || for the solution of the 
equation of motion of inverted pendulum. Such an ansatz (7)-(8) is physically justified, because due to inertia the 
system responds weakly to fast external pulses. 

Substituting (7), (8) into (1), (2) and expanding all functions in power series of p and retaining terms up to the 

first order (e.g. gf } (q, p) ~ gf '(Q, P) + p( Xj • V Qj )gf } (Q, P) + p( Pj ■ V Pj )gf ) (Q, P), etc.), we obtain 

dH . _ s dH , _ .dH gl (0) (Q,P) 

M% • V Qj ) g -° )( ^ P) + • Vp 3 )^M + £ g f ) (Q , P) cos(fcm) + 

k>l 

A*E^ • V Qj )gf ) (Q,P)cos(fcm) +/i^( /0j Vp 3 )gf ) (Q,P)cos(fcm); (9) 

fc>i fc>i 

dH . „ . di? , „ .dH hf } (Q,P) 
Pi +„ Pi = - /*(* ■ Vq 3 .)— - ■ Vp 3 )— + + 

M*i • V Qj ) h '' 0(Q ' P) +M(ft • V Pj ) h ' (0) ( 2 Q ' P) + ^ h| fc > (Q, P) cos(fcm) + 

fc>i 

A* • V Q, ) h f° (Q. p ) cos(fc^) + /* J>,- V P , )h< k) (Q, P) cos(fcfit). (10) 
fc>i fc>i 

(Here we use the Einstein summation convention on repeated indices). All functions in these equations depend on 
mean variables (Q, P). The equations contain slow and fast terms which can be separated by averaging over the period 
t = 2ir/n. More precisely, (9) and (10) are of the form: Slow + Fast = 0. Since fl ^> uiq, the slow terms may be 
considered constant over the time 2tt/Q and therefore may be substituted by their averages. The equations take the 
form: (Slow}n + Fast ~ 0. On the other hand, it is evident that {Fast)n = and consequently the result is a set of 
coupled equations for the slow and fast variables: {Slow)n — 0; Fast = 0. Hence, the first step consists in integrating 
(9), (10) with respect to time over the period 27r/f2. Notice that Q, P, as well as H(Q, P, t) are considered to remain 
constant when integrating over "fast time ". Only slow terms withstand this operation. The remaining terms form 
the fast equations. Particular attention is to be given to the last and next to the last terms of (9), (10). These terms 
contain the products X j cos(kflt) and Pj cos(fcili). Inasmuch as X j an d pj may contain both cos(fc'ili) and sin(fc'i7i) 
(as will become apparent later), we have here the products cos(k'ilt) cos(kttt) and sin(fc'ftt) cos(fcfM). The average of 
the first product is nonzero for k' — k, whereas the average of the second product always vanishes. In other words, 
the last two terms in (9), (10) may in general contribute both to slow and fast equations. These latter thus read 
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The terms in (11), (12) do not all have the same order. While the terms /x X~ ~ u oX > M P~ M^P ~ w oP are not 
small, the terms with /jx, MP are much smaller. So, we find for the fast equations accurate up to zero-order terms in 
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Integrating these equations, we obtain 
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Subsequently, we solve the fast equations taking into account the terms of first-order in zx too, substituting (14) into 
(11), (12) and again integrating with respect to time (recall that Qi, Pi, H(Q,P,t) are held constant) 
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Expressions (15), (16) represent a rapidly oscillating response of coordinates and momenta of the system to the 
fast external actions. The point is that it is not an unique reaction of the system. The mean trajectories (Q, P) are 
altered too. Inserting (15), (16) into (9), (10) and averaging over the period 2n/Q, we obtain the following equations 
of motion for the mean canonical coordinates: 
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It is pertinent to note that all additional terms (~ l/^ 2 ) in these equations come from the last two terms of (9), (10) 

27r/f2 

when integrating of the function cos(fc'f^) cos(kQt) over time: J cos(fc'Oi) cos(kflt)dt — hSkk' ■ 

o 

Relations (17), (18) constitute the effective equations of motion of a system subject to high-frequency external 
actions. The respective Liouville equation for the effective distribution function F(Q, P,t) follows immediately from 
(17), (18) in the form (4)-(5): 
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where Q = (Q, P) is the phase space vector of the mean coordinates, J is the inverse Jacobian of the transformation 
from the initial mean phase coordinates to the coordinates at time t: 
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The phase space velocity vector Q is given by the right-hand sides of the equations (17), (18). 

In the simple case of a forced nonlinear oscillator with a Hamiltonian H = p 2 /2m + U(q) and the driving force 
.F(g)cosfit, we have the only nonvanishing amplitude h^(q) = F(q) and from the equations (17), (18) follows the 
equation of motion in the mean coordinate 
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This coincides with the result first obtained by Kapitza Q, || . The last term in (21) stabilizes the inverted pendulum. 
In other cases, there may be different extra terms in equations (17), (18) and respectively in the Liouville equation. 
These terms depend on various combinations of derivatives of the functions g- fe \ h- fe ^ with respect to canonical variables 
and may be responsible for many interesting dynamical and thermodynamic effects. 

It is important to emphasize that even in the case when the equations (1), (2) are Hamiltonian, the equations 

(17), (18) may be non-Hamiltonian, i.e. the compressibility does not in general vanish x(Q) — Vg 0. It is not 
surprising since a driven system behaves like a thermodynamically open system. 

We would like to thank Dr. M.E. Tuckerman for helpful correspondence on the subject matter of this paper. This 
work was supported by Conselho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq, Brazil). 
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